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Abstract 

Theoretical studies on M13 (M = Fe, Co, Ni) and Mi3Pt n (for n = 3, 4, 5, 20) clusters including 
the spin-orbit coupling are done using density functional theory. The magnetic anisotropy energy 
(MAE) along with the spin and orbital moments are calculated for M43 icosahedral clusters. The 
angle-dependent energy differences are modelled using an extended classical Heisenberg model with 
local anisotropies. From our studies, the MAE for Jahn- Teller distorted Fei3, Mackay distorted 
Fei3 and nearly undistorted C013 clusters are found to be 322, 60 and 5 /xeV/atom, respectively, 
and are large relative to the corresponding bulk values, (which are 1.4 and 1.3 //eV/atom for 
bcc Fe and fee Co, respectively.) However, for N113 (which practically does not show relaxation 
tendencies), the calculated value of MAE is found to be 0.64 //eV/atom, which is approximately 
four times smaller compared to the bulk fee Ni (2.7 /ieV/atom). In addition, MAE of the capped 
cluster (Fei3Pt4) is enhanced compared to the uncapped Jahn- Teller distorted Fei3 cluster. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



Density-functional theory (DFT) is an adequate tool to obtain useful information on the 
physical properties of small clusters. However, with respect to magnetic properties, the 
variations of magnetic moments with cluster size and morphology are often too obscure for 
establishing a clear trend (for example, see the behavior of the exchange coupling of free Fe 
clusters in [Hl2])- Still, certain general statements can be made. It may be safely concluded 
that the local magnetic moments in outer shells of clusters are enhanced, as compared to the 
interior of clusters or to the corresponding bulk crystal [3], IH El E] . This is a typical effect 
of reduced atomic coordination at the surface [7] , well known from numerous experiments 
or calculations of magnetic slabs. Moreover, an analysis of cluster morphologies, by both 
experiment and theory [HI |9|, [TUJ [TTJ [TJl [13j [HI [15] , reveals the importance of icosahedral-like 
structures, prohibited for bulk or surface phases. On the other hand, the diversity of cluster 
structures in combination with the surface enhancement of magnetic moments make clusters 
interesting model objects for tuning magnetic properties at the nanometer scale. 

When addressing magnetic properties of clusters from the computational point of view at 
a realistic level, it is important to take into account two specific issues: (i) a possible non- 
collinear (NC) setting of magnetic moments (i.e., a smooth variation of magnetic density 
vector from point to point in space), and (ii) spin-orbit interaction (SOI), along with the 
existence of orbital moments. Both issues have a long record of incorporation into first- 
principles DFT calculations, and are internally related: They mix the spin-up and spin-down 
states and must be, in principle, treated alongside on equal footing (see, e.g., Refs. jTHJ [TTj for 
a review). Discussing specifically DFT calculations for clusters, one notes certain technical 
difficulties in combining NC spin density with the SOI, which acted so far as a limiting 
factor on the number of calculations performed, and the size of clusters treated: Lack of 
symmetry, big effect of structure relaxation, slow convergence, large size of simulation cell 
around a free cluster. Moreover there is a conceptual problem of choosing a "correct" 
non-collinear solution among many apparently close metastable configurations. Different 
groups report very different results for the same systems (e.g., Refs. [TBI EES EQ]), so that a 
preference of one or another result is not obvious. Therefore, not so much the size of cluster 
is a problem in itself, as the organization of a calculation in a way allowing to extract 
meaningful results, that is, clearly defining structural and magnetic models and carefully 
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analyzing their consequences. 

A better understanding of the origins and details of large orbital magnetic moment as 
well as large MAE in clusters [2IJ [22] is demanding in order to manipulate materials for 
improving the technology like in magnetic data storage devices. Binary 3d— 5d clusters can 
be a challenging material in this respect. For example, there are observations that Ad and 
5d elements like Rh, Pt and Au, which are nonmagnetic in bulk phase, attain significant 
moments when alloyed with 3d transition metals like Fe, Co or Ni [23l |24"1 25, 26J. 

In the present work, we analyze the effect of capping icosahedral clusters M 13 of M = 
Fe, Co, and Ni with Pt atoms. Theoretically, there have been several studies on magnetic 
anisotropy for supported clusters [271 ESI ESI ESI EH [32], [33], whereas studies on small free 
clusters are still limited [311 EHl [361 E3 [38]. The motivation for the present study is that 
icosahedral clusters for Mi 3 are known to be very stable, [HI [121 ES] and the alloying of 3d 
transition metals with Pt results in large magnetic anisotropy as well as orbital moment. 
Hence, it is expected that the orbital and the spin moments of M 13 clusters be strongly 
affected by capping with Pt. 

Most of our DFT calculations have been done with the planewave code Vasp (Vienna 
Ab-initio Simulation Package) [101 [HJ H2]. For test purposes and for validating the results 
of calculations on non-capped icosahedral clusters, we performed calculations also with a 
local-orbitals code, Siesta (Spanish Initiative for Electronic Simulations with Thousands 
of Atoms) [43] . As the two methods are very different in what regards the technical imple- 
mentation of the DFT calculation scheme, and both have been earlier used in calculations 
of magnetic clusters, their direct comparison might be of its own interest. In addition, few 
test calculations for the binary clusters were done using the all-electron local-orbital code 
Fplo (Full-Potential Local-Orbital scheme) [44J. In this respect, it is interesting to note 
that all-electron calculations confirm the results obtained with Vasp and Siesta. 

The paper is organized as follows. Section [IT] outlines the calculational methods 
and setup. Section III deals with the results for monometallic icosahedral clusters, 
notably comparison between Vasp and Siesta results. Section IV discusses the re- 
sults for capped clusters, obtained with the Vasp code. We have made a comparison 
between Vasp and Fplo for one of the capped clusters. Conclusions are drawn in Section V. 
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II. COMPUTATIONAL METHODS 



Most of the calculations have been performed with the Vasp code [lUl HH E2] based 
on DFT and within the generalized gradient approximation (GGA). The parameters by 
Perdew are used for the exchange and correlation functional |45]. Vasp uses the projector 
augmented wave (PAW) method (42j |46] and a planewave basis set. Periodic boundary 
conditions were imposed onto large enough cubic cells (with the edge of 15 A for M13 
clusters and 20 A for M 13 Pt ra clusters) in order to sufficiently reduce the interaction between 
replicated cluster images. Only the T point was used for the Brillouin-zone sampling for the 
cluster calculations. A k-mesh of (llxllxll) is used for the bulk calculations to compute 
the equilibrium lattice constants in case of bcc Fe, fee Ni and Co. The values for local 
magnetic moments correspond to the integration of continuous magnetization density over 
atom-centred spheres with radii of 1.302 A (Fe, Co), 1.286 A (Ni) and 1.455 A (Pt). It 
is a well known fact that the magnetic anisotropy energy for cubic bulk transition metals, 
defined as the maximum energy difference, per atom, between different settings of the spin 
moment relative to the atoms framework, is of the order of 10~ 6 eV. Moreover, a special care 
is demanded for the study of the MAE in clusters because any slight error can accumulate 
and produce misleading results while dealing with energy differences of such a small scale. 
In order to calculate the magnetic anisotropy of M13 icosahedral clusters, we have used the 
energy convergence criterion for the self-consistency as 10~ 10 eV with a Gaussian half-width 
parameter of 0.01 eV for the discrete energy levels. A very high plane wave cutoff value of 
1000 eV as well as a dense Fourier grid spacing of 0.04 A in each of x, y and z directions is 
taken. 

While relaxing the clusters using conjugate gradient algorithm, the same energy conver- 
gence criterion is used but a lower plane wave cut-off of 270 eV. The structural relaxation of 
clusters are done in the scalar relativistic mode. For the calculation of orbital moment and 
MAE, the SOI is treated as a perturbation in the scalar relativistic Hamiltonian. Bench- 
mark calculations for some M13 clusters have also been done with the Siesta code [43] 
within GGA, where the exchange and correlation functional is parameterized by Perdew, 
Burke and Ernzerhof |4Tj . It uses the norm- conserving pseudopotentials of Troullier and 
Martins [48J. The localized basis functions of "double-^ with polarization orbitals" quality 
(and triple-^ for 3d functions) have been constructed according to the standard scheme of 
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the Siesta method [15] , with the "Energy Shift" parameter of 0.01 Ry. The treatment of 
SOI was included as described by Fernandez- Seivane et al. [30]. In the Siesta calculations, 
we used the atomic coordinates as relaxed by Vasp. This allows for a direct comparison 
of non-collinear structures, spin and orbital moments. The standard representation of local 
magnetic moments in Siesta is in terms of Mulliken populations of localized orbitals, that 
is a quite different definition from that used in Vasp. For the sake of better comparison, 
we report in the following the magnetic moments as properties integrated over atom-centred 
spheres for both Vasp and Siesta. 

TABLE I: Relaxed coordinates (A) for Jahn- Teller distorted (JT) and Mackay transformed (MT) 
Fei3 clusters. See Fig. 1 for an illustration and labelling of atoms. 







Fei 3 (JT) 






Fei 3 (MT) 




Atom Nr. 


X 


y 


z 


X 


y 


z 


1 


0.6764 


-2.0818 


-1.0306 


0.0000 


-1.2917 


-2.0157 


2 


-0.6764 


2.0818 


1.0306 


0.0000 


1.2917 


2.0157 


3 


0.0000 


0.0000 


-2.3361 


0.0000 


1.2916 


-2.0157 


1 


0.0000 


0.0000 


2.3361 


0.0000 


-1.2916 


2.0157 


5 


-0.6764 


-2.0818 


1.0306 


-1.2917 


-2.0157 


0.0000 


6 


1.7709 


-1.2866 


1.0306 


1.2916 


-2.0157 


0.0000 


7 


-1.7709 


1.2866 


-1.0306 


-1.2916 


2.0157 


0.0000 


8 


0.6764 


2.0818 


-1.0306 


1.2917 


2.0157 


0.0000 


9 


-1.7709 


-1.2866 


-1.0306 


-2.0157 


0.0000 


-1.2917 


10 


2.1890 


0.0000 


-1.0306 


2.0157 


0.0000 


-1.2916 


11 


-2.1890 


0.0000 


1.0306 


-2.0157 


0.0000 


1.2916 


12 


1.7709 


1.2866 


1.0306 


2.0157 


0.0000 


1.2917 


13 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 



Finally, several test calculations have been performed by a third method in order to 
confirm the presented results. These calculations are documented later in Section IV and 
have been carried out with the all-electron local-orbital code (Fplo 6.00-24) [H] in its clus- 
ter mode [51] using LSDA [52]. The valence basis comprised 3s3p3d 4s4p4d 5s states for 
3d-metals, while 5s5p5d 6s6p6d 7s states were used for Pt. The fully relativistic mode is 
employed here, where Fplo solves the four-component Dirac-Kohn-Sham equations includ- 
ing spin-orbit coupling (at the single electron level) in all orders [53J. Default settings were 
used for the other numerical parameters. 
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III. RESULTS FOR MONOMETALLIC ICOSAHEDRAL CLUSTERS M 



In order to systematically pursue a search towards probable relaxation patterns from 
the ideal icosahedron (ICO) structure, we "drove" the structure along two different paths, 
which are known from previous studies [121 EH ESI EE] and are likely to lead to different 
metastable arrangements. The structural relaxation of M13 clusters results in Jahn- Teller 
(JT) distortion as well as partially Mackay transformed (MT) clusters (Fe, Co, Ni), which 
is shown in Fig. 1 for the case of Fei 3 . For the Fei 3 cluster, the JT-distorted structure is by 
123 meV/cluster, the MT cluster by 35 meV/cluster lower in energy compared to the energy 
of the ideal ICO, see also Ref. [H]. Only the Fe cluster exhibits the JT-distortion and the 
Mackay distortion as additional local minimum in the energy curve [MJ EE] , whereas Co and 
Ni show only a tendency for MT. 

What we refer to in the following as the Jahn- Teller (JT) type distortion [58] is that 
which maintains the fivefold symmetry around one of the ICO axes; it allows a compression 
or tension of the cluster along this axis, possibly accompanied by a mutual opposite rotation 
of two pentagonal rings pierced by the axis in question. 

The occurrence of the JT distortion originates from (accidental) quasi-degeneracy of 
the highest occupied molecular orbitals so that the system may lower its energy by level 
splitting arising from the induced distortion and corresponding redistribution of electrons. 
The partially Mackay distortion is of completely different origin since it is connected with 
the tendency of the magnetic Fe clusters to transform to the fee cuboctahedron (CUBO) 
and by subsequent Bain like transformation to bec Fe. Whereas for large Fe clusters, the 
JT distortion is no longer observed, we find in the simulations that the partially Mackay 
transformed cluster (which were denoted as shell- wise Mackay transformed cluster in [X2|l56]) 
still corresponds to a metastable [59] (local energy minimum) state up to cluster sizes having 
15 closed atomic shells (n) and magic atom numbers (N) defined by 

N = |(10n 3 + 15n 2 + lln + 3). 

For the Fei3 cluster, the JT distortion corresponds to a compression of the cluster along 
the depicted in Fig. 1, whereas we have chosen a different orientation for the 

MT cluster in the same figure (z-axis through mid-points of two opposite bonds) in order 
to illustrate its more cubic-like appearance. The relaxed atomic coordinates for the MT 
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TABLE II: The bond lengths in A for relaxed Mackay transformed M13 clusters. 



Bonds 


Fei 3 


Co 13 


Nil 3 


Ccntrc-shcll 


2.394 


2.33 


2.321 


24 X 


2.501 


2.451 


2.441 


6 X 


2.583 


2.447 


2.439 



and JT 13-atom Fe clusters are listed in Table I. For the JT-distorted Fei3 cluster, 10 of 
the peripheral atoms lie at a distance of 2.42 A from centre, the other two atoms show an 
inward relaxation towards the centre and hence, remain at a distance of 2.34 A from centre. 
A clear picture of the relaxation is shown in the left panel of Fig. 1, where the arrows 
associated with atoms on the top and bottom of the box face oppositely. The structural 
relaxation of the Fei3 cluster shown in the right panel in Fig. 1 corresponding to the partial 
Mackay transformation is the onset of the transformation from an ICO to CUBO, where 
the consecutive triangular facets of the ICO change to square facets of the CUBO, which 
can be described by a parameter s defined as the square of the ratio of the stretched to 
the unstretched edges being 1 and 2 for the ICO and CUBO, respectively. See [57] for the 
original explanation. It is well observed for multishell ICO-like Fe clusters [T2l 155] . In our 
calculations, the relaxed Fei3 has s = 1.07; C013 and M13 have s ~ 1, resulting in an almost 
ideal icosahedral like structure. 

The peripheral bond lengths for the MT clusters (see Table II) can be arranged in groups: 
6 of the 30 bonds being parallel to the cartesian axes (x: 5-6, 7-8; y: 1-3, 4-2; z: 9-11, 10-12; 
see Fig. 1 for the labeling of the atoms), the remaining 24 bond lengths are also listed in 
Table II. Remarkably, the distances of the 12 peripheral atoms to the central one are all 
identical. 

The relaxed MT coordinates obtained with Vasp have been used in the calculation by 
Siesta, without further relaxation, in order to compare the resulting values of spin and 
orbital moments, and the non-collinearity. Table III compares the calculated results by the 
two methods for Fei3(MT) and C013 (MT) clusters. For both methods, the spin moments 
remain parallel (x- and y-components of the spin vectors are, at most, 0.001 //b, and are 
not shown in Table III). The orbital moments for MT distorted Fei3 obtained by both 
programs are directed towards the z— direction in the lower part of the cluster and are 
deviated "outwards" from the z— direction in the upper part, while for C013, they show 
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TABLE III: Spin and orbital moments of the relaxed Fei3 (MT) and C013 clusters as calculated by 
Vasp and Siesta, for the initial [001] setting of spin direction. The x and y-components of the spin 
moment are < 0.001, hence not shown here. For Siesta calculations, the values of spin moments 
are given in terms of Mulliken populations, the magnetic moments are in units of /is/atom. 



Vasp Siesta 



Fei 3 


Atom No. 




hy 


L z 


S z 


La; 


hy 


L z 


s z 


1 


0.00 


-0.02 


0.11 


3.06 


0.00 


-0.01 


0.09 


3.43 


2 


0.00 


0.02 


0.12 


3.06 


0.00 


0.01 


0.09 


3.43 


3 


-0.01 


0.00 


0.08 


3.08 


-0.01 


0.00 


0.08 


3.43 


4 


0.01 


0.00 


0.08 


3.08 


0.01 


0.00 


0.08 


3.44 


5 


0.00 


0.00 


0.07 


3.10 


0.00 


0.00 


0.06 


3.44 


6 


0.00 


0.00 


0.07 


3.10 


0.00 


0.00 


0.06 


3.44 


7 


0.00 


0.00 


0.07 


3.10 


0.00 


0.00 


0.06 


3.44 


8 


0.00 


0.00 


0.07 


3.10 


0.00 


0.00 


0.06 


3.44 


9 


0.01 


0.00 


0.08 


3.08 


0.01 


0.00 


0.08 


3.44 


10 


-0.01 


0.00 


0.08 


3.08 


-0.01 


0.00 


0.08 


3.44 


11 


0.00 


0.02 


0.11 


3.05 


0.00 


0.01 


0.09 


3.43 


12 


0.00 


-0.02 


0.11 


3.06 


0.00 


-0.01 


0.09 


3.43 


13 


0.00 


0.00 


0.05 


2.70 


0.00 


0.00 


0.03 


2.75 


C013 


1 


0.00 


0.03 


0.09 


2.10 


0.00 


0.01 


0.08 


2.42 


2 


0.00 


-0.03 


0.09 


2.10 


0.00 


-0.00 


0.08 


2.41 


3 


0.02 


0.00 


0.14 


2.10 


0.02 


0.00 


0.10 


2.42 


4 


-0.02 


0.00 


0.14 


2.10 


-0.02 


0.00 


0.10 


2.42 


5 


0.00 


0.00 


0.17 


2.10 


-0.00 


-0.00 


0.13 


2.41 


6 


0.00 


0.00 


0.17 


2.10 


0.00 


-0.00 


0.13 


2.42 


7 


0.00 


0.00 


0.17 


2.10 


0.00 


-0.00 


0.13 


2.42 


8 


0.00 


0.00 


0.17 


2.10 


0.00 


-0.00 


0.13 


2.41 


9 


-0.02 


0.00 


0.14 


2.10 


-0.02 


0.00 


0.10 


2.42 


10 


0.02 


0.00 


0.14 


2.10 


0.02 


0.00 


0.10 


2.42 


11 


0.00 


-0.03 


0.09 


2.10 


0.00 


-0.00 


0.08 


2.41 


12 


0.00 


0.03 


0.09 


2.10 


0.00 


0.01 


0.08 


2.42 


13 


0.00 


0.00 


0.05 


1.91 


0.00 


0.00 


0.04 


1.96 


opposite 


behavior, i.e. 


the orbital moments deviate 


away from the 


z— axis in 


the lower part 



of the cluster and they are directed towards the z— axis in the upper part, but we note 



that these deviations are not significantly above the numerical noise level. This becomes 
evident if one compares data for those atoms that are mutually equivalent (apart from 
rotations around the z-axis) by symmetry: 1-4, 5-8, 9-12. While the symmetry requirement 
is (almost) obeyed by the atoms 5-8, deviations of a few 1/100 fi B are found between the 
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TABLE IV: Cluster averaged values for (L) = ^£i=i |Lj| and (S) = ^£i=il s il of MT M 13 
clusters are shown in /is/atom, compared to bulk values for bcc Fe, fee Co and Ni as obtained 
from our calculations. 



Cluster (L) 


<s> 


(ILDbulk 


(|S|> bulk 


Fei 3 0.084 


3.04 


0.06 


2.25 


C013 0.125 


2.08 


0.08 


1.67 


Nii 3 0.070 


0.66 


0.05 


0.65 


TABLE V: The values of prefactors in the Heisenberg 


; anisotropy term 


in Eq. (2) for the relaxed 


M13 clusters. 








Cluster 


D 2 (meV) 


D 4 (meV) 


D 6 (MeV) 


Fe 13 (ICO) (46 Mb) 






-41.1 


Fei 3 (JT) (44 Mb) 


15.2 






Fe 13 (MT) (44 p B ) 




-11.5 


-41.1 


C013 (31 mb) 




-21.7 


-25.1 


Nii 3 (8 mb) 




7.1 


7.3 



atoms 1-4 and also between the atoms 9-12. 

When comparing the numerical results for spin and orbital magnetic moments from these 
two different calculational methods, one must take into acount the difference in their defini- 
tion. In Vasp, the properties (spin and orbital moments) are extracted as projection onto 
an atomic sphere (see Section [Tl]). The "standard" Vasp value of atomic sphere radius for 
Co does in fact correspond to slightly overlapping spheres in our C013 cluster. The Siesta 
output results are reported in terms of decomposition over projection onto localized, but 
spatially extended, numerical orbitals, known as Mulliken population analysis. It is known 
that the local magnetic moments as well as atomic charges in heterogeneous systems do 
often come out very different, when estimated according to these two different schemes. In 
order to illustrate this effect, we give in the last column of Table III the values of spin mo- 
ment, extracted from the Siesta results. The fluctuations of these integrated values over 
apparently equivalent atoms are caused by the sparseness of the spatial grid with the step 
0.078 A, on which the summation of spin density has been done. 

The (L) and (S) of Fei 3 , C013 and Ni 13 clusters are compared with respect to bulk in 
Table IV. The orbital and spin moments for bulk systems are calculated with the equlibrium 
lattice constants of 2.832 A (bcc Fe), 3.52 A (fee Co) and 3.523 A (fee Ni). The orbital and 
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spin magnetic moments of the elemental Fe, Co and Ni clusters are larger than the related 
bulk values. For example, in Ref [60J it was shown that in Ni clusters of up to 13 atoms, 
the average orbital moment (L) per atom is found to be 4 to 8 times larger with respect 
to Lbuik; for larger cluster sizes, (L) was shown to approach the bulk value. It turns out 
from our calculation that for NI13 clusters, (L) agrees quite well with the earlier report [60J, 
with orbital moment having larger magnitude than the bulk value. For Fei3 and C013, the 
(L) is also relatively large with respect to the bulk systems. In addition, the average spin 
moments (S) for these clusters are also increased with respect to the bulk systems. 

In addition, the magnetic anisotropy energy for M 13 clusters was calculated. Bulk 3d 
metals show a tiny (cubic Fe, Co and Ni) MAE of ~1 /xeV. The tiny value in cubic systems 
is due to the high symmetry and it is expected that an ideal ICO exhibits a similarly small 
MAE. If the symmetry is broken, e.g., by tetragonal distortion of a cubic system, the MAE 
increases considerably (see for instance Ref. [61]) . 

Figure 2 shows the ideal icosahedral cluster and the definition of the (x, z)-plane with 
angle 9 for the magnetization directions used in our MAE calculations. The magnetic 
anisotropy energy is computed from the total energy as AE = E tot (8) — E tot (6 = 0), E tot (8) 
being the ^-dependent total energy. The ^-dependent energy obtained from total energy 
calculations is compared with the anisotropy term of the extended Heisenberg model, defined 
as 

H — — JijSiSj + H ani + //dipolar, (1) 

where 

TV 

H ani = -J2J2 D ^- S ^ n - ( 2 ) 

n i 

The spin-spin and dipolar interactions are neglected, while comparing ab-initio total 
energy differences with the Heisenberg model. 

In Eq. (2), n is the order of anisotropy (an even integer), are the unit vectors of 
the atoms along the radial directions, Sj their spin moments and D is the anisotropy 
constant (of Neel-like model), which is assumed to be ^-independent. The second and 
fourth order contributions are constant from symmetry considerations of a ideal ICO. The 
anisotropy consists then mostly of a sixth order contribution, assuming that still higher 
orders are negligibly small. Figure 3 shows a comparison of energy differences of an ideal 
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Fei3 ICO obtained from ab-initio calculations to that of the Heisenberg model using a 
least-mean square fit to the ab-initio data. It is a delicate task to correctly determine 
energy contributions in the order of a few /ieV/atom or below. To obtain meaningful results 
presupposes that the charge density is extremely well converged and still consistent with 
the symmetry of the system. A distortion of the symmetry might occur due to numerical 
fluctuations (which can to a limited extent be controlled by, e.g., the choice of energy cutoff 
and Fourier-grid) and to perturbations by the setup of the model system itself, e.g., due 
to the electrostatic interaction between periodic images of the supercell. Suitable settings 
can be found from the corresponding contributions to the forces, which depend on the ionic 
positions and the charge distribution as well. Consequently, the magnitude of the forces 
acting on all equivalent atoms should be nearly equal, which in our case, we have taken 
care that forces are not larger than 10~ 6 eV/A. It is shown from model calculations of fee 
Co clusters that the ^-dependence of magnetic anisotropy may also differ with the polar 
angle [H2], which is not discussed here. 

Along with the MAE for a ideal cluster shown in Fig. 3, the MAE of relaxed M13 clusters 
is calculated: The variation of energy difference with respect to 9 for relaxed clusters are 
shown in Fig. 4. The ^-dependent energy difference for the JT-distorted Fei3 cluster (the top 
panel of Fig. 4) gives rise to a significant D 2 contribution of ~ 15.2 meV for the anisotropy. 
However, approximately only a D 4 and D 6 contribution arises for the anisotropy in case 
of MT Fei3 cluster. The values are depicted in Table V. The MAE for the JT-distorted 
Fe cluster is calculated to be 322 /ieV/atom, which is approximately five times larger with 
respect to the MT cluster (MAEm^ = 60 /ieV/atom). The reason behind the large MAE 
for the JT-distorted cluster relative to the MT cluster may be attributed to the breaking of 
symmetry in the former case. 

The ^-dependent energy difference for the relaxed MT Fei3 cluster shows qualitatively 
good resemblance with that of the ideal cluster (see Fig. 3), however, the quantitative MAE 
value for ideal Fei3 ICO is found to be approximately 1 /xeV/atom, while for the MT Fei3 
cluster, it is approximately 60 times larger than for the ideal cluster. However, for relaxed 
geometries of C013 and Ni 13 , the MAE value is calculated to be 5 and 0.64 /ieV/atom, 
respectively. Our calculations suggest a larger value of MAE per atom for Fei 3 and C013 
compared to the bulk, which is 1.4 /ieV/atom (bec Fe), 1.3 /xeV/atom (fee Co) [S3]- For Ni 13 , 
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FIG. 1: The JT-distorted (left) and the MT (right) Fei3 cluster. Arrows indicate the direction 
of relative shift of atoms with respect to the ideal positions. For the JT-distorted and MT Fei3 
cluster, the displacements marked by the arrows have been scaled up by a factor of 20 and 30, 
respectively. The box is only guide to the eyes. (The actual simulation box size is 15 3 A 3 .) 

this is lower with respect to bulk fee Ni (2.7 //eV/atom). The small anisotropy for the relaxed 
Nii3 cluster relative to bulk Ni is due to the high symmetry of the cluster, M13 having a 
nearly ideal icosahedral shape. The change of sign in the energy vs. ^-curves (shown in Fig. 
4) for the relaxed M13 clusters can be explained through their way of structural relaxation. 
Namely, the MT relaxation pattern for the C013 cluster is just opposite to that of Fei 3 : 
The Co "dimers" oriented along the edges of the cube in the right panel of Fig. 1 make the 
shortest side in each C03 face, whereas the corresponding Fe-Fe bonds are the longest in 
the Fe3 faces of the MT Fei3 cluster, see Table II. By comparing ab-initio results with the 
Heisenberg anisotropy term, the D4 of M13 can be fitted, which are listed in Table V. It is 
observed that D4 makes a significant contribution to the magnetic anisotropy of the relaxed 
clusters having cubic-like distortions. 

It would be fascinating to see how spin and orbital moment as well as the MAE change, 
if we consider more "asymmetric" clusters. To this we have added a varying number of Pt 
atoms on the top of M i3 clusters. 
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FIG. 2: The ideal icosahedral cluster showing the x-z plane in which the angle 9 is varied in the 
MAE calculations for the partially Mackay transformed clusters. 

IV. RESULTS FOR CAPPED CLUSTERS Mi 3 Pt n 

Obviously, binary and core-shell transition metal clusters show a yet larger diversity 
compared to the elemental systems. For example, the intermixing of 3d elements with 
4d or 5d elements results in a large magnetic moment of the binary systems [26]. Both 
in experiment and in calculations, it has been observed that FePt and Co clusters show 
enhanced spin moments and orbital moments with respect to corresponding bulk values [641 
163] . Hence, it would be interesting to study the change in magnetic properties including 
MAE of M13 clusters capping with Pt atoms, which is described in the following. 

We considered three high-symmetry positions to cap the M13 clusters by a single Pt atom 
and found that a Pt position above the centre of a facet is most favorable in all three cases. 
In the following, we used this finding as a guideline for initial geometries of Mi3Pt n clusters 
(n = 3, 4, 5, 20): In all cases, the Pt atoms were initially placed above the facet centres at 
a distance found in the single-Pt capping case. After relaxation, optimised geometries were 
obtained as shown in Fig. 5 for Coi3Pt3 (left) and Coi3Pts (right) and in Fig. 6 (right side) 
for Nii 3 Pt2o- For n = 20, the initial geometries form a core-shell morphology. 

Since atom projected quantities like spin and orbital moments depend on the specific 
code, we compared for the particular case of Nii3Pt3 related data obtained by Vasp and 
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FIG. 3: The ^-dependent energy for an ideal Fei3 icosahedral cluster showing the ab-initio data 
(filled circles and solid line as guide to the eyes) compared with the Heisenberg model (dashed 
curve). For Vasp calculations, the energy difference is in units of meV per cluster and is defined 
as AE = E — E e=0 - 

Fplo. The structure optimisation was carried out by VASP and the same geometry was 
used to evaluate the magnetic moments by both codes. 

In Fplo, the magnetic moments are calculated through Mulliken population analysis. 
Figure 7 shows the absolute value of orbital moment per atom |Lj| (left) and the absolute 
value of spin moment per atom |Sj| (right) as a function of distance of each atom from 
the cluster centre for Ni 13 Pt 3 . It is obvious that both codes give results which are in good 
agreement with each other. Therefore, the calculations which are discussed in the following 
are done with Vasp only. 

In Figs. 8 and 9, the variation of |L$| and of |S$| with respect to the distance from centre of 
Mi 3 Pt n clusters are shown, respectively (the symbols are kept consistent for both figures and 
the centre atom is placed at zero). Figure 8 nicely shows that the orbital moments of such 
few-atom systems in general depend very sensitively on the particular chemical composition 
and geometry. The resulting electronic structure can be very individual (e.g., the nature 
of the highest occupied level depends on the electron number and the spin moment) and is 
hard or impossible to predict without a detailed calculation. It is also found that (see Fig. 
8) | Li | for M13 clusters approach towards the corresponding bulk values of bcc Fe, fee Co 
and Ni reported in Table IV, with increasing the number of Pt atoms. The orbital moment 
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FIG. 4: The plots from top to bottom show the energy differences in meV/cluster vs. 8 of relaxed 
Fei3 (f° r two types of relaxations - JT and Mackay distortion), C013 and M13 clusters. In case 
of JT-distorted Fei3: 9 varies with respect to the z-axis in the plane passing through atom 4, 
middle of the bond between atoms 6-12, atom 10 and atom 3 (shown in the left side of Fig. 1), 
whereas for the MT cluster, 9 varies through all atoms as depicted in Fig. 1. Energy difference is 
(AE = E — Eg=o). The energy difference for the C013 cluster is multiplied by a factor 10 and the 
energy difference for the N113 cluster is multiplied by a factor 100. 

on the centre Fe atom for Fei 3 Pt n cluster increases with the number of Pt atoms, whereas 
for Co 13 Pt„ and Ni 13 Pt n , this trend is weak. The orbital moment on Pt atoms is found to 
be very sensitive with respect to the core atomic species. For example, the |Lj| on Pt atoms 
for Nii 3 Pt n clusters are enhanced relative to that of Fei 3 Pt„ and Coi 3 Pt n . 

In Fig. 9 the variation of spin moments for the capped clusters are shown, which suggests 
that with increasing Pt atoms, the |Sj| for M 13 clusters always remain larger with respect 
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FIG. 5: The relaxed structure of Coi3Pt3 (left) and CoisPts (right) shows the arrangement of Pt 
atoms on the C013 cluster, dark (blue) spheres: Co atoms and light (yellow) spheres: Pt atoms. In 
each case all Pt atoms lie in the same plane corresponding to the lowest energy structures found 
so far. 




FIG. 6: The initial and final strucutres for the Nii3Pt2o cluster are shown in the left and right panel, 
respectively. Ni and Pt atoms are represented by dark (blue) and light (yellow) spheres, respec- 
tively. The final structure, corresponding to the lowest energy structure found so far, demonstrates 
the importance of atomic relaxations. 

to the corresponding bulk 3d-spin moments values. Unlike the trend in orbital moment, the 
spin moment on Pt atoms is not much affected by its core species. In the following, a few 
quantitative statements are made for each of the capped clusters. 
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FIG. 7: The variation of the individual atomic |L*| (left panel) and |Sj| (right panel) with respect 
to the distance from centre for Nii3Pt3 cluster. : results by Vasp and Q: results by Fplo. The 
open and filled symbols (for both colors) show the moments on Ni and Pt atoms, respectively. 

A. Fei3Pt n clusters 



The values of average orbital moment and average spin moment for each atomic species 
along with the total orbital moment L tot and total spin moment Stot for the capped clusters 
are defined and reported in Table VI. Our observation suggests that adding more Pt atoms 
on Fei3 cluster, the L to t for Fei3Pt n and S tot for Fei3Pt3 and Fei3Pt5 increase considerably. 
However, there is a decrease of S tot for Fei 3 Pt 2 o because of the antiferromagnetic like align- 
ment of spin of centre Fe atom with respect to the sorrounding atoms in this cluster (for 
all other cappings, the Fei3 cluster is ferromagnetic). Capping with three and five Pt atoms 
enhances the atomic orbital moment on Fei3 by ~ 50 % and 33 %, respectively, relative to 
the bulk value of bcc Fe, i.e. 0.06 /is/atom. It does not much change the Fe moment in 
comparison with the bare Fei3 cluster, however. The average orbital moment (Lm) on Fei3 
decreases monotonically with increasing number of Pt atoms. 
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FIG. 8: The variation of the individual atomic orbital moments with respect to distance from 
centre for Fei3Pt n (left), Coi3Pt n (middle) and Nii3Pt n (right) clusters. The light (green) and 
dark (black) symbols for all cases represent the |Lj| on Pt atoms and M43 clusters, respectively. 

B. Coi3Pt„ clusters 



A ferromagnetic ground state is found to be stable for all compositions of Coi3Pt n clusters, 
which leads to a monotonically increasing trend of L tot and S to t 111 these clusters. The (Lm) 
and (Lpt) for Coi3Pt n clusters shows a trend similar to Fei3Pt n , i.e. with increasing number 
of Pt atoms, the (Lm) (as defined in Table VI) on the C013 cluster decreases and the (Lp t ) 
on the Pt atoms is enhanced. The (Sm) on C013 is merely constant, while it increases 
monotonically for Pt with more Pt atoms on the cluster surface. 
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FIG. 9: The plots in the left, middle and right panels show the variation of |Sj| with respect to the 
distance from centre for Fei3Pt n , Coi3Pt n and Nii3Pt n clusters, respectively. The same symbols 
are used here as in Fig. 8. 

C. Nii3Pt„ clusters 



In this case, the L tot and S to t also show a regular increase with the number of Pt atoms. 
The (Lm) and (Lp t ) for this cluster do not much depend on the number of Pt atoms. 
However, a decrease of (Lp t ) and (Spt) from Nii 3 Pt5 to Nii 3 Pt2o cluster is also observed. 
This is probably caused by a structural instability occuring for this composition upon the 
relaxation. The geometry optimization of this cluster converges to a structure with different 
symmetry, with a ferromagnetic ordering, where the Ni atoms are placed closer to the surface 
of the cluster as shown in Fig. 6. The reason for the segregation of Ni atoms towards the 
surface may be due to its lower surface energy compared to Pt [66] . Another aspect related 
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TABLE VI: The orbital and spin moments for the binary Pt-M clusters in ^s/atom, where we 
have distinguished core and shell contributions. Cluster averaged values are defined for (Lm) = 
]^Si=i and (Sm) = jsYli=i |Sz| (M represents the 13-atom Fe, Co and Ni clusters). The 
averaged values for (Lp t ) = ^J27=i 1-^*1 an d (Spt) = ^Y17=i ( n ^ s * ne num ber of Pt atoms and 
here, n= 3, 5, 20). |L to t| and |Stot| are the corresponding absolute values of the total orbital and 
spin moment for every composition of the binary clusters. 



Cluster 


(Lm) 


(Lpt) 


<S M > 


(Spt) 


|Ltot| 


|Stot| 


Fei 3 Pt 3 


0.09 


0.08 


2.86 


0.36 


1.43 


38.3 


Fei 3 Pt 5 


0.08 


0.134 


2.97 


0.526 


1.76 


41.3 


Fei 3 Pt 2 o 


0.05 


0.17 


2.18 


0.401 


2.74 


36.4 


Coi 3 Pt 3 


0.1 


0.13 


1.87 


0.42 


1.68 


25.6 


Coi 3 Pt 5 


0.095 


0.16 


1.83 


0.47 


2.06 


26.1 


Coi 3 Pt 2 o 


0.07 


0.3 


1.88 


0.60 


4.9 


36.2 


Nii 3 Pt 3 


0.06 


0.2 


0.81 


0.402 


1.61 


11.7 


Nii 3 Pt 5 


0.05 


0.26 


0.80 


0.512 


2.01 


13.0 


Nii 3 Pt 2 o 


0.08 


0.2 


0.81 


0.36 


3.55 


17.6 



to this may be observed in the right panel of Figs. 8 and 9 showing the variation of onsite 
orbital and spin moments with respect to the distance. The large variations in orbital and 
spin moments just occur because of the structural distortion for this cluster composition. 
Comparing all three cases of capped clusters, we find that the presence of Pt atoms on M13 
affects the orientation of core orbital moments in such a way that they always prefer to be in 
non-collinear alignment for the M 13 Pt n clusters, which is not the case in the uncapped M 13 
clusters. On the other hand, directions of individual spin moments remain always collinear 
for the same clusters indicating that they are less affected by the Pt atoms. 

Regarding the MAE of capped Mi 3 Pt n clusters, we observe that the symmetry rules the 
magnitude of the effect like for M13 systems. For example, in Fei3Pt4 (see Fig. 10) with 
second-order anisotropy, we find a high MAE value, even larger by a factor of two (~ 7 
meV/cluster) compared to JT-distorted Fei 3 . 
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9 (degree) 

FIG. 10: Left: The x-z plane (where 9 is varied) of the Fei3Pt4 cluster, showing the positions of 
the Pt-atoms relative to the ICO in the unrelaxed configuration. Light (orange) and dark (blue) 
spheres are the Fe and Pt atoms, respectively. Right: The ^-dependent energy differences for the 
relaxed Fei3Pt4 cluster. The energy difference is in units of meV per cluster and is defined as 
AE = E — Eq = q. 

V. CONCLUSION 

The spin and orbital magnetic moments are calculated for M13 and Mi3Pt n clusters 
including spin-orbit interaction using DFT methods. The studies on MAE of M 13 clusters 
show that the calculated MAE values of relaxed M 13 clusters are considerably enhanced 
relative to both the ideal icosahedral clusters and to the corresponding bulk values in case 
of Fei3 and C013. For Nii3, the MAE is smaller than that of the bulk value of fee Ni. The 
MAEs for relaxed clusters are found to be affected by the degree of relaxation due to which 
Fei3 clusters have a high anisotropy in comparison with C013 and Nii3. Iron is a special 
case for which we obtain two local energy minima [T2l [56] . corresponding to the JT and MT 
cases. The JT-distorted Fei3 cluster exhibits an anisotropy approximately 5 times larger 
than the partially MT cluster due to the low symmetry in the former case. The present 
calculations of MAE agree well with the qualitative prediction of the Heisenberg model for 
the ^-dependent energy differences. With respect to the spin and orbital moments, both the 
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capped and the free clusters show an increased value of orbital and spin moments compared 
to the bulk. For the capped clusters, the spin moments on Pt atoms remain unaffected by 
the host atoms. Finally, we infer that deposited transition metal clusters, with very large 
effects of relaxations, may even exhibit larger MAE values. Self-assembly of such clusters, 
like in case of Fe-Pt [13] and Co [H], may then approach the class of functional magnetic 
materials of use for magnetic storage devices. 
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